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We study conformational and electrophoretic properties of polyelectrolytes (PEs) in tetrava- 
lent salt solutions under the action of electric fields by means of molecular dynamics simulations. 



C4H ' Chain conformations are found to have a sensitive dependence on salt concentration Cs- As Cg is 

o ■ 

■ increased, the chains first shrink to a globular structure and subsequently reexpand above a critical 



concentration C*. An external electric field can further alter the chain conformation. If the field 
strength E is larger than a critical value E* , the chains are elongated. E* is shown to be a function 



Q ' of Cs by using two estimators Ej and Ejj through the study of the polarization energy and the 

o : 

onset point of chain unfolding, respectively. The electrophoretic mobility of the chains depends 

1^ , strongly on Cs , and the magnitude increases significantly, accompanying the chain unfolding, when 

' E > Ejj. We study the condensed ion distributions modified by electric fields and discuss the 

O ■ 

l/-^ • connection of the modification with the change of chain morphology and mobility. Finally, E* is 

, studied by varying the chain length N. The inflection point is used as a third estimator Ejjj. Ejjj 

scales as A^^'^-^3(^) and A^^0 ''6(2) at Cg = 0.0 and C*, respectively. Ejj follows a similar scaling law 
to Ejjj but a crossover appears at Cs = C* when N is small. The Ej estimator fails to predict the 
critical field, which is due to oversimplifying the critical polarization energy to the thermal energy. 
Our results provide valuable information to understand the electrokinetics of PE solutions at the 
molecular level and could be helpful in micro/nano-fluidics applications. 
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I. INTRODUCTION 



Electric field-driven molecular analysis and sorting techniques have been widely used in 
many domains of research, such as chemistry, biology, and medical engineering. While the 
functionality of a bio-sensing or diagnostic system becomes more and more complicated and 
specialized today, electrokinetics remains the mechanism of choice for fluid actuation and 
manipulation at micrometer or submicrometer scales through the use of electric fields 
Therefore, a solid understanding of the electrokinetic behavior of ions, molecules, and macro- 
molecules under the influence of electric fields is necessary to successfully integrate electric 
fields in micro/nano-fluidic devices. 

Electrophoresis has been developed to separate charged macromolecules, such as DNA 
molecules or proteins, based upon their molecular weight for many years {2 4|. In most 
situations, electrophoresis is performed in sieving mediums such as gels. This is because, in 
free solutions, DNA molecules cannot be size-separated owing to the free-draining effect: the 
hydrodynamic friction and the molecular charge of a DNA molecule are both linearly pro- 
portional to the chain length. Hence, the electrophoretic mobility is length independent jsl. 
This phenomenon applies not only to DNA molecules but to any charged macromolecule, or 
polyelectrolyte (PE). In gels, the mechanism of electrophoretic separation is dominated by 
biased reptation of PEs and so the mobility has chain-length dependence {3,6]. One draw- 
back of the method is its low efficiency. Researchers continue to search for new techniques 
to separate charged molecules more efficiently j^, particularly in micro/nano-fluidics. 

Recently, Netz proposed an idea to separate PEs electrophoretically in free solutions [?, 8|. 
PEs are polarized by electric fields and, if the electric field is strong enough, the polarization 
can induce chain unfolding, which renders a drastic increase in the electrophoretic mobility 
of the chains. Most importantly, the critical electric field to unfold chains is found to depend 
on the length. This dependence provides a plausible manner to separate chains by size in free 
solutions through an unfolding transition. If the idea could be applied in micro/nano-fluidics, 
the efficiency to separate DNA or other biomacromolecules by size could be improved. In 
Netz' idea-demonstrating work j?, 8|, PEs were simply collapsed by monovalent counterions 
by setting a strong Coulomb coupling parameter between charged particles. However, this 
situation is not realistic because the PEs are usually collapsed by adding condensing agents 
such as multivalent salt {q]. Coulomb coupling is not as strong as in the simulation work 
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and furthermore, there are several species of counterions present in solutions which compete 
with each other to condense on the chains. The distributions of these ions play a crucial role 
in determining how PEs unfold in electric fields. Simulating polarization induced unfolding 
in salty solutions requires that individual ions be explicitly modeled with interactions that 
realistically account for aqueous conditions. This is what we will focus on in this study. 
When salt is added in solutions, especially multivalent salt, PEs show complicated be- 



havior 
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3Nl2|. The addition of multivalent salt can induce collapse or aggregation of chains, 
which causes phase separation. An excessive addition of multivalent salt can even cause the 
separated phases to dissolve back to a homogeneous solution; the collapsed and aggregated 
chains reexpand and separate from each other. Using these phenomena, researchers are able 



to control the dimension o 
ordered toroidal particles 



DNA molecules in solutions and collapse them into small, very 



13|. In addition to the size control, the presence of multivalent 

n 

counterions can also modify the charge distribution around PE chains [1^ and lead to a 
specific phenomenon called "overcharging", where the condensed counterions overcompen- 



sate the charge on chain surfaces 11 



Iq . The total chain charge is also effectively altered. 



In certain conditions, the effective chain charge changes sign, a phenomenon called "charge 



inversion ' 



12| . A simple way to determine the effective charge is to study electrophoresis of 



16|. Controlling the charge of PE-ion complexes is an important 



PEs in weak electric fields 
issue for gene therapy because it is related to the efficiency of DNA up-take by cells through 
endocytosis pathways 1^]. 

In strong electric fields, the behavior of PEs becomes even more complicated because the 
charged particles in the complexes respond to the field in different ways, depending on the 
charge and the sign. Moreover, strong electric fields can change the molecular conformation 
in ways that are difficult to predict. The electrokinetics of the system is significantly mod- 
ified. It has been demonstrated that DNA molecules condensed by polyvalent counterions 
such as spermine can be decondensed in DC electric fields if the electric field strength ex- 



ceeds a threshold E* 
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Netz predicted that E* should scale as '^'^^'^ where N is the 



chain length and z/ is the Flory exponent [7|, Isj. Simulations did observe scaling behavior 



but the scaling exponent measured did not agree with the prediction 



19|-|2l|. Hsiao and Wu 



considered a coiled chain as an ellipsoidal object of volume V and proposed a modification 
of the sc aling law to be E* ~ V~^^'^ 19|. The salt valency dependence of E* has also been 



explored 



19 



21| . The results show that the magnitude of electrophoretic mobility of chain 



3 



increases significantly above a chain unfolding transition, providing the foundation for chain 
separation in free solutions. Moreover, the chains can be unfolded to an elongated structure, 
which can then be utilized in the techniques of single-DNA molecule sequencing to increase 



spatial resolution of detection 22 



The response of PEs to alternating-current (AC) electric fields has also been investigated 



recently by simulations 2l|,|25|. Liu et al. found that chains are stretched and the sizes breath 



with the frequency of applied AC field only when the field strength exceeds some critical value 
and the frequency is smaller than the intrinsic relaxation frequency of the chain 2l|. The 
work by Hsiao et al. further connected the critical AC field strength with the DC one, and the 



critical AC frequency with the inverse DC chain- fluctuation time 25|. A model, based upon 
Maxwell- Wagner dielectric theory, has also been developed, which explains the critical field- 



frequency correlation for chain unfolding in AC fields [25|. Recently, the conformational 
transition to a stretched state in AC fields has been experimentally demonstrated 26 |. 
The chain size shows interesting hysteretic behavior upon sweeping the AC frequency. To 
understand the behavior, it is very important to first investigate chain conformations, ion 
distributions, and also the mobilities of the chains and ions in DC fields, because DC fields 
can be regarded as AC fields with zero frequency. 

PE solutions involve both polymeric and electrolyte degrees of freedom, which brings 
many difficulties in dealing with these systems theoretically. Molecular dynamics simu- 
lations are a simple and economic tool, able to study systems in a controllable way and 
capture detailed information at the molecular level. Many simulation works investigate the 

[ioi,Q,0jr 



structure of PEs in different solution conditions 
fields 0, S 
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30| and the response to electric 
32| . However, detailed information about the electrokinetics 



of different species of ions around the chains is lacking. Since the conformation of chains 
depends strongly on the concentration of salt in solutions, it is very important to understand 
the distribution of ions around the chains and see how these ions are affected by external 
electric fields. Therefore, in this study we aim to investigate electrophoresis of single PEs in 
multivalent salt solutions by means of molecular dynamics simulations. The rest of the paper 
is organized as follows. We explain the model and simulation setup in Sec. II. The results 
and discussions are presented in Sec. HI. The effect of tetravalent salt concentration on chain 
conformation is discussed first (Sec. HI- A). We then study polarization and determine the 
critical electric field at various salt concentrations (Sec. III-B). The electrophoretic mobility 
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of chains and condensed ions are investigated in Sec. III-C. The distribution of condensed 
ions and the effective charge of chains are presented in Sec. III-D and Sec. III-E. Finally, 
the mobility dependence on chain length is studied in Sec. III-F. We give our conclusions in 
Sec. IV. 



II. MODEL AND SIMULATION METHOD 



Our system comprises a single linear chain, modeled by a bead-spring chain model. The 
chain consists of monomer beads. Each bead carries a negative unit charge — e and 
dissociates one monovalent cation (or "counterion" ) into the solution. The bonds connecting 
two adjacent monomers are modeled by the finitely extensible nonlinear elastic (FENE) 
potential 



f/FENE(fe) = -2<a.ln [I - j j (1) 

where b is the bond length, b^^x is the maximum extension, and k is the spring constant. 
Salt is added into the system. The salt molecules dissociate into tetravalent cations (also 
called "counterions" ) and monovalent anions ( "colons" ) in the solution. All the particles — 
including the monomers, counterions, and colons — are modeled explicitly as spheres that 
are described by the purely repulsive Lennard- Jones (LJ) potential 

ULj{r) = \ - (2) 

1^ , for r > ^(j 

where r is the separation distance; a and ^lj represent the diameter and the hardness of 
the LJ sphere, respectively. Since the interaction between monomers is purely repulsive, our 
system corresponds to a good solvent. Charged particles also interact with each other via 
the Coulomb interaction 

f/c(r) = kBT\B^ (3) 
r 

where T is the temperature, is the Boltzmann constant, and Zi is the charge valency of the 
ith particle. = e^/ {AireQerkBT) is the Bjerrum length where eo is the vacuum permittivity 
and is the relative dielectric constant of the solvent. At the separation distance r = Ag, the 
electrostatic energy between two unit charges is exactly the thermal energy ksT . The solvent 
molecules are not modeled explicitly in the study. However, their effects are incorporated 
implicitly through the following three ways: (1) the dielectric constant e^, which takes into 



account of the dielectric screening of charge in the solvent medium, (2) the friction force 
—rriiQvi, which models the drag acting on particle i, proportional to the moving velocity, 
(3) the stochastic force 57j(t), which simulates the thermal collisions of solvent molecules on 
the particle i. The equation of motion is the Langevin equation, 

niiri = -miCiVi + + ZicEx + rfj (4) 

where mj is the particle mass, Q is the friction coefficient, and = —dU/dfi is the 
conservative force. In Langevin dynamics, the temperature is determined by the fluctuation- 
dissipation theorem: {f]i{t) ■ ffjit')) = QkBTmiQ5ij5{t — t'). The external electric fleld is 
uniform and exerts an electric force ZicEx on the particle i in the x-direction The system is 
placed in a periodic rectangular box. Particle-particle particle-mesh Ewald sum is used to 



calculate the Coulomb interaction 
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We assume that all the particles have the identical mass m and LJ parameters a and ^lj- 
We set Slj = 0.8333/cbT, k = 5.8333A;BT/cr^, bma.x = 2(T, = 3a, and Q = lr~^, where r = 



aWm/(kBT) is the time unit. The chain length is varied and the monomer concentration 



Cm is flxed at 0.0003cr~^, which describes a dilute polymer solution 3^. The simulation 
box is adjusted linearly with in x-direction, and has a volume of 1.3A^ x 50.64 x 50.64. 
This prevents chains from self-overlapping through the periodic boundary condition. The 



salt concentration Cs is varied from = to C^_ = 0.0006cr ^, a range for the chain to 



exhibit the behavior of reentrant condensation 



28 



430l]. The electric fleld strength E is varied 



over a wide range, from E = to 2kBT/ (ea), to study its effect on the properties of the 
PE system. The Langevin equation is integrated by the Verlet algorithm. The integration 



time step At is equal to 0.005r 35|. A pre-run phase takes about 10'' time steps to bring 
the system into a stationary state, followed by a production-run phase of 10® time steps to 
cumulate data for analysis. Since h ydro dy namic interaction is largely screened out under a 



typical electrophoretic condition we neglect the hydrodynamic interaction m this 

study. Recent simulations have demonstrated the validity of this approximation when the 
chain length is not short (A^ > 20) [3^. We study the electrophoretic properties of chains 
at a flxed chain length A^ = 48 for the flrst 5 subsections of Sec. Ill whereas in Sec. III.F, A^ 
is varied from 12 to 768. To simplify the notation, the value of a physical quantity will be 
reported in the {a, m, ksT, e)-unit system in the rest of the text. For example, the strength 
of electric fleld is described in units of ksT/^ea), the dipole moment is in units of ea, and 
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the electrophoretic mobility is in units of ea'^ /{rkBT), and so on. 



III. RESULTS AND DISCUSSIONS 



A. Chain conformation in electric fields 



In this section, we study the conformation of PEs at different salt concentrations Cg under 
the action of electric fields. The chain length is 48. We first calculated the mean square 
radius of gyration R^, which is used to characterize the chain size. The definition is given 
by = Xlili (i^i " ^cmy) /N where is the position vector of the monomer i and fcm is 
the center of mass of the chain. The results are plotted in Fig. [1] Each curve shows how R^ 
varies with Cs at a given field strength E. 
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FIG. 1: Mean square radius of gyration as a function of tetravalent salt concentration Cs at 
different strength E of electric field. The value of E can be read in the figure. The error bar of 
data in this paper is smaller than the size of data symbol if it is not presented. 



We see that in zero electric field, R^ decreases with Cs up to the salt concentration 
C* = 7.5 X 10~^. The value of C* is one fourth of Cm, at which the amount of tetravalent 
cations in the solution are in charge equivalence with the monomers on the chain. The 
decrease of R^ shows that the chain collapses, which is due to ionic screening and to the 
bridging effect induced by the added tetravalent salt. At Cs = C*, the chain is nearly 
neutralized by the condensation of tetravalent counterions. In the region Cs > C*, R^ 
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becomes an increasing function of Cg- The chain size reexpands shghtly. The collapse and 
the reexpansion of the chains can be regarded as a single-chain version of the PE precipitation 
and redissolution. These phenomena are collectively called "reentrant condensation" 15|. 
Reentrant condensation for PEs has been investigated in detail in theories 
experiments 
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15 



40 



43 



44l ]. and simulations |27H30|. 



Now let us focus on the cases in which the electric field is applied. All these curves 
lie above the zero- field limit. This result demonstrates that the chain is elongated by the 
applied electric field and therefore, takes a larger value of Rg. The stronger the electric field, 
the larger the deviation from the zero-field limit. However, this deviation also depends on 
Cg- While increased with E at a given Cg in the region Cg < C*, is basically unchanged 
with E in the region Cg > C* for E < 0.4. When E > 0.5, the PEs no longer exhibits a 
compact structure at C*. The chain size continues to shrink with Cg for Cg > C*. This 
shrinkage suggests that the excess of the tetravalent counterions in the solution helps the 
chains to collapse, against the action of the electric fields. In a very strong field such as 
E = 2.0, R^ becomes basically a constant with Cg. The chain collapse by tetravalent salts 
is completely suppressed by the strong electric field. 

We next calculated the asphericity A of the chain in electric fields. The 
quantity describes the degree of geometrical deformation away from a sphere and 
can be used to characterize the chain conformation. It is defined hj A = 
i(((Ai-A2)2 + (A2-A3)2 + (A3-Ai)2)/(Ai + A2 + A3)^) where Ai, A2, and A3 are the 
three eigenvalues of the gyration tensor of the chain. The tensor was calculated by 
Tai3 = ~ rcm}a{ri — rcm)ii/N, whcrc the subscripts a and /3 denote one of the 

three Cartesian components x, y, z of the subscribed vector, respectively. The value of A 
ranges from to 1. It is for a perfect sphere and 1 for a rod. For a random coiled chain, 
A takes the value 0.431 obtained by simulations 45|. The left panel of Fig. [2] presents the 
asphericity A of our chain as a function of the electric field E. Each curve denotes one case 
running at a given Cg. 

The curves are horizontal lines when E is small, showing that chain distortion does not 
happen in the weak electric fields. At Cg = 0.0, the value of A is 0.68. It is larger than 
0.431, showing that the chain expands more than a coil does. We found that A decreases 
with Cg in this weak E'-region, and reaches a minimal value 0.20 at Cg = C*. The small 
value of A shows the formation of a compact globule structure. When Cg > C*, A begins 
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FIG. 2: Asphericity A and degree of chain unfolding Ds vs. electric field strength E at different 
Cs- The value of Cs is indicated in the legend. The symbol "0" in the right panel denotes the 
threshold electrical field (denoted by Ejj in the next section) at a given Cg- 

to increase with Cg. The result demonstrates a chain reexpansion. In the intermediate E 
region between 0.1 and 1.0, A exhibits a drastic increase, showing the deformation of the 
chain by the electric field. At large field E = 2.0, the value of A is around 1. Therefore, the 
chain is deformed from its "natural" conformation in the zero field to a rodlike structure. 

In the right panel of Fig. [21 we show the degree of chain unfolding Dg as a function of E 
at different Cg- Here Dg is defined to be the ratio of the chain end-to-end distance Re to the 
chain contour length L^. Similar to A, Dg displays a sigmoidal increase when electric field 
is above some threshold value. In very high electric fields, the value of Dg can be as large 
as 0.9, showing that the chain is quasi-fully stretched. We estimate the threshold field by 
the onset point, at which Dg increases 10% from its zero-field limit. The dependence of the 
threshold field on Cg can be clearly seen in the figure. 

To support the results of our calculation, we present in Fig. [3l Panel (a), (b), and (c), 
the snapshots of simulation at three salt concentrations Cg = 0, 7.5 x 10~^, and 6.0 x 10~^, 
respectively. The three Cg present the three cases with the amount of the adding salt smaller 
than, equal to, and larger than the equivalence point C*. 

The left picture in each panel of the figure shows the case subjected to a weak electric 
field E = 0.02, which is below the critical field E*. We can see that the chain exhibits an 
expanded-coil structure in the salt-free solution (Panel (a)), a compact globule structure at 
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(a) 





(c) 

FIG. 3: (Color online) Snapshots of simulation at (a) Cs = 0.0, (b) Cg = 7.5 x 10~^, and (c) 
Cs = 6.0 X 10"^. In each panel, two pictures are presented: the left one is the case in a weak 
field E = 0.02 and the right one is in a strong field E = 1.0. The PE is presented by a yellow 
bead-spring chain. The monovalent counterions, the tetravalent counterions, and the monovalent 
colons are presented by green-colored, red-colored, and white-colored spheres. The direction of 
electric field points toward the right direction. 
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the equivalence point (Panel (b)), and a less compact structure at the high salt concentration 
(Panel (c)), consistent with the calculations of the asphericity in Fig. |2j The right picture in 
each panel presents the case in the strong electric field E = 1.0, which is above E*. We can 
see that the chain unfolds to an elongated structure. The chain size in Panel (a) is longer 
than in Panel (b), and than in Panel (c), which agrees with the results in Fig. [T] where 
Rg is a decreasing function over Cg for E = 1.0. Moreover, we can see that the elongated 
chain is aligned parallel to the field direction and the condensed monovalent counterions 
(in Panel (a)) and the condensed tetravalent counterions (in Panels (b) and (c)) cumulate 
more densely near the chain end which follows the field direction (the right-hand end in 
Fig. 13] figure). It shows that the chain is polarized in the electric field and possesses a dipole 
moment pointing to the field direction to reduce the system energy. 

This dipole moment is precisely what causes the chain to unfold. The electric field 
causes an inhomogeneous distribution of counterions to condense on the chain. Due to this 
inhomogeneity, the electric force acting on the two sides of the chain is not balanced, which 
gives an effective tension along the chain. The chain is unfolded if this tension is strong 
enough to break the electrostatic binding between monomers by the condensed counterions. 

B. Polarization and critical electric field 

The previous section discussed the existence of a critical field E*, beyond which a PE 
is aligned along the field direction and drastically unfolded. To understand the unfolding 
mechanism in details, consider the polarization of a chain of polymerization = 48 in 
electric fields. 

The energy of the electric field stored in a dielectric material is = (1/2) J D-EdV where 
D is the electric displacement field. Since D = e^E + P where P is the polarization density, 
the energy W can be written as a sum of two terms: the first is Wq = (1/2) J eo\E\'^dV, which 
represents the field energy in vacuum, and the second is Wpoi = (1/2) j P ■ E, which denotes 
the polarization energy. In our study, the polarization energy is calculated by Wpoi = p- E/2 
where p is the induced dipole moment of the PE-ion complex. For an unfolding event 
to occur, the polarization energy Wpoi is larger than the thermal fiuctuation energy ksT, 
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according to Netz 7|, |8|. The dipole moment of the PE complex can be calculated by 



P= ^ Zie{fi-rcm] 

iSPE— complex 



(5) 



where fi is the position vector of the zth particle and i runs over the entire PE complex 
(both the monomers and the condensed ions), and rem is the center of mass of the entire 
complex. There exists no definite way to define a PE complex. In this study, we primitively 
define the PE complex by a constant threshold distance r^: a PE complex comprises the 
chain itself and the ions with the distance to the chain smaller than r^. These ions are said 
to condense on the chain. We chose = A^- This criterion has been used previously to 



study ion condensation on a rigid chain 46| but other definitions can be used 39|. Since 
the polarization will occur in the field direction, we calculated here the a:-component of the 
dipole moment, Px- The results are presented in Fig. HJas a function of E for different Cg. 
The criterion for chain unfolding is that > ^ksT/E, and a dashed line demarcates the 
boundary of the two regions in the figure. Wpoi is larger than ksT above the line, whereas 
smaller below the line. 
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FIG. 4: Dipole moment as a function of E. Each curve corresponds to a salt concentration 
whose value can be read in the figure. The dashed line is the equation px = 2kBT/E* and the 
dotted hne is a reference hne with the slope equal to 1. 



The points of Px fall on a straight line in the log-log plot in the weak field region E < 
2kBT/px. The slope of the line is 1, which follows exactly linear response theory, px = aE. 
The polarizability a is directly related to the height of the line in the figure, a takes a large 
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value for a non-collapsed chain, and is smallest for a collapsed chain at Cs = C*. In the field 
region above the dashed line, px deviates from this simple power law. Approximately above 
the dashed line, the curves grow faster than the linear dependence but then curve down, 
exhibiting a hook-like behavior. 

In the weak field region, the chain conformation is not perturbed by the applied electric 
field. However, the displacement of the condensed ions on the chain is still possible. It is 
hence the ion displacement that is responsible for the formation of the dipole moment and 
results in the linear response behavior. When the electric field is strong enough to unfold the 
chain, the chain elongation provides the further possibility for the condensed ions to migrate 
on the chain 25|. Consequently, the dipole moment acquires a value larger than predicted 
by linear response theory. In very strong electric fields, the chain reaches its maximum 
extension. Further elongation of the chain becomes impossible. Increasing the electric field 
can cause condensed ions on the chain to be stripped off. The number of the condensed ions 
decreases. As a consequence, the dipole moment decreases, resulting in the hook-like curve. 

The critical electric field E* is estimated as the intersection of the simulated px curves 
with the unfolding transition px = 2kBT/E*. The obtained critical field is reported in Fig. [5] 
as a function of Cs and denoted Ej (for reasons that will be clear momentarily). We can see 
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FIG. 5: Critical electric field E* as a function of Cs- E* is determined by two different methods: 
£^1 from the polarization energy and E'Jj from the onset point of Dg. 



that Ej increases with Cg and reaches the maximum value at Cs = C*, and then, decreases 
slightly. This behavior is consistent with what we have observed in Fig. [H The chain is 
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easier to stretch in the low-salt region {Cs < C*) than in the high-salt region {Cg > C*). 
To verify that the chain begins unfolding aX E = Ej, the threshold field obtained from the 
onset increasing of Dg in Fig. |2]is plotted and denoted by Ejj. The consistency between Ej 
and Ejj data seems to suggest that chain unfolding occurs as when linear response theory no 
longer holds. Nonetheless, further verification through varying the chain length from = 12 
to 768 shows that Ej is not always situated at the break-down of the linear dependence, as 
shown in Fig. [6l 




FIG. 6: Dipole moment px as a function of -E at Cs = and Cg = C*. Chain length is indicated 
in the right panel. The dashed line defines the equation px = 2kBT/E* . The symbol "0" denotes 
the dipole moment at E = Ejj. 



Indicating the onset-critical point on the curve, we clearly see that the dipole moment 
deviates from the linear behavior at E = Ejj, but not at E = Ej. Therefore, the criterion 
Wpoi > ksT for chain unfolding is too simple to be precise. A real unfolding point can 
happen with the polarization energy larger or smaller than ksT, depending on both the 
chain length and salt concentration. It necessitates a fundamental understanding of the 
unfolding mechanism in the future to set up a correct criterion for the problem. The results 
shown here demonstrate that chain unfolding is tightly connected with the polarization 
change. 
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C. Electrophoretic mobility 



We have seen that PE-ion complexes exhibit a drastic unfolding transition when the 
electric field is stronger than E*. One pertinent question is whether the mobility of the chain 
shows a drastic change too, accompanying the conformational transition. Also, it is very 
important to know the mobility of ions, especially the condensed multivalent counterions, 
because these ions play a decisive role in determination of the chain conformations. 

The mobility of the chain and of the condensed tetravalent counterions are defined as 
fJ'pe = Vpe/E and fi']^^^ = v\^/E, respectively, where Vpe and are, in turn, the velocities 
of the chain and the condensed tetravalent counterions in the field direction. The results 
are presented in Fig. O Panels (a) and (b), for /ipe and /i!^4, respectively, as a function of E 
for different Cs- The sign of the obtained mobility can be positive or negative, denoting the 
moving direction of the studied object toward the +x- or the — x-direction. To investigate 
the relationship between the mobility change and the conformational transition, we have 
plotted on the curves of the data the corresponding critical field E* (the E*jj) by the symbol 
"0". 

We can see in Fig. [Tl^a) that when E < E*, constant and depends on the salt 

concentration. For Cs < C*, the value of fipe is smaller than zero because the effective chain 
charge is negative, the same sign as the bare chain charge. Increasing Cg leads to a decrease 
of \fipe\- The absolute net chain charge decreases due to the condensation of the counterions 
on the chain backbone. At Cs = C*, fipe is approximately zero. This is because the total 
charge of the tetravalent counterions in the solution is equivalent to the charge of the PE. 
The chain is effectively neutralized upon the condensation of these ions. If Cs is increased 
above C*, fipe becomes positive. The chain now moves in the +x-direction and the effective 
charge of the PE-ion complex is positive. In other words, a charge inversion occurs. Charge 



inversion induced by multivalent salt has been observed in experiments 



47|. Nonetheless, 



researchers continue to put their efforts on this study area for a full understanding of the 



underlying mechanism [16|, |48|, |49[. When the electric field is strong E > E*, fipe is no 



longer a constant but monotonically decreases with E. The fact that the effective chain 
charge becomes more negative suggests that a strong electric field strips the condensed 
counterions off the chain. A detailed study concerning the number of the ions condensing on 
the chain will be presented in the next section. Our results reveal that charge inversion can 
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FIG. 7: Mobility of (a) chain fipe and (b) the condensed tetravalent counterions function 
of E. Each curve denotes one case running at a specific Cs- The value of Cs can be read in the 
figure. The critical field Ejj is indicated on the corresponding curve by the symbol "0". (c) fipe 
and ^5^4 vs. Cg in different E fields, replotted from the data of (a) and (b). The value of E is 
indicated in the left panel. 



be suppressed by strong electric fields. Comparing mobility with the conformational change 
demonstrates that chain mobility and unfolding are closely related. 

Fig. El^b) shows the mean mobility of the tetravalent counterions fi'^^ condensed on the 
chain. We observed that /i^^4 is approximately equal to fipe when E is small. We recall that 
the mobility of tetravalent counterions is positive in a free solution, but Fig. Wlh) shows 
that /i!j_4 < when Cg < C*. A negative /i!j_4 clearly demonstrates that the electrostatic 
interaction between tetravalent counterions and monomers sets a strong constraint on the 
ions. The ions are dragged along with the chain and consequently, /i^4 ~ fipe. When E 



16 



is strong enough to break the constraint, jj,'^^ deviates from /ipe- In very strong electric 
fields, these ions glide along the chain surface in the +x-direction; fi'^^ thus takes a positive 
value. The mobility difference between the chains and the condensed ions suggests that ion 
condensation must take place in a dynamic way. The condensed ions glide on the chain 
and are eventually stripped off the chain end. The ions from the bulk solution continually 
condense onto the chain from the other end. The process repeats continuously and a balance 
is established. 

Please notice that the field required to cause /i^^4 to deviate from /ipe is stronger than 
E*. This suggests that the chain unfolding occurs before the ions are able to glide along the 
chain. Therefore, there is a field region in which the chain has unfolded but the condensed 
tetravalent counterions are too tightly bound to glide. Since fipe becomes more and more 
negative, so does /i^j_4 before shooting up in stronger E fields. Consequently, the /i^4 curve 
exhibits a minimum. For Cs > C*, the binding between condensed tetravalent counterions 
and monomers becomes effectively weaker because of overcharging (see Fig. [8]^b)). The field 
region where the chain unfolds without ion gliding shrinks as Cg is increases. The minimum 
disappears and /i^j_4 grows monotonically with E. Fig. [7](c) shows the mobility presented as 
a function of salt concentration in different field strength. As we can see, fipe and yU^4 both 
increase with Cs- Moreover, /Xpe is roughly equal to /i^j_4 over the studied Cs for E < 0.4. For 
E > 0.4, very different behavior is observed. The mobility of the tetravalent ions increases 
much faster than /ipe, and eventually becomes entirely positive at very strong E while /ipe 
remains negative. 

We remark that in this study, each data point of the mobility is calculated under the ac- 
tion of an electric field at a given strength, once the system reaches a steady state. Figs. [Tl^a) 
and (b) are not obtained by sweeping the DC field strength at a constant rate. Therefore, no 
hysteresis is displayed. Our results report the mobilities in a stationary, pseudo-equilibrium 
condition. For simulations done by sweeping the field strength, one would expect the occur- 
rence of a hysteretic behavior in mobility, and also in chain size, if the sweeping rate is fast. 
There exists two relevant characteristic times in chain conformational transitions: the first 
is the chain unfolding time to transition from a coiled or globular chain to an elongated one, 
and the other is the collapse time to relax from an elongated chain to a coil or globule. The 
unfolding time is generally longer than the collapse time because additional time is needed 
to disentangle a coiled or globular chain. The chain size (and therefore mobility) will follow 
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a different path when sweeping from a weak DC field to a strong one and then from a strong 
field to a weak one, presumed that the total sweeping time would be comparable to the 
two characteristic times. This deserves further investigation in the future. A counterpart 
has been shown in experiments studying chain-size hysteresis when AC field frequency is 
changed at a fixed field amplitude {26!] . 



D. Number of condensed ions and effective chain charge in electric fields 

In this section, we study the number of ions condensed on the chain {N = 48) and effective 
chain charge in electric fields. As in the previous section, an ion is regarded to condense onto 
the chain if its distance to the chain is smaller than the threshold distance rt = Xb- Once 
the condensed ions are identified, the effective chain charge can be calculated by summing 
the charges of the condensed ions and the chain monomers. Since there are three types of 
ions in the simulation box, we treat each species independently. The results are presented 
in Fig. [H] as a function of E where A^^^, N^^, and N^-^ are the numbers of the condensed 
monovalent counterions, tetravalent counterions, and colons, respectively. In order to make 
comparison with the chain conformational transition, we have plotted the chain asphericity 
A at three representative Cg in the figure. 

We observe that N'i^^ is unperturbed by E, when E < E*. This constant value decreases 
with Cs and goes to zero when Cg > C*: There are no condensed monovalent counterions 
above the equivalence point. When E > E*, N'^^ shows a decreasing behavior with E for 
Cs < 3.75 X 10"^ but at = 5.62 x 10"^ and 7.5 x 10-^ the A^^^ exhibits a small hump. For 
even higher Cs, the curve increases slightly. By comparing with the chain asphericity, we 
find that these variations take place at the moment when the chain changes its conformation. 
Therefore, the number of condensed ions is closely related to chain morphology. How salt- 
induced chain conformation affects the ion condensation in zero electric field has been studied 
in reference 
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N^i at Cs = 0.0 can be used to verify Manning condensation theory [46j. The theory 
states that counterion condensation takes place on a rigid PE while the mean distance per 
unit charge on the chain is smaller than the Bjerrum length. It results in an effective chain 
charge equal to —cLc/Xb- In our simulations, the mean distance before condensation takes 
place is equal to the averaged bond length (b) = 1.1. Hence Manning's theory predicts an 
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FIG. 8: Number of condensed (a) monovalent counterions N^-^_, (b) tetravalent counterions 
and (c) coions N^i as a function of E. The salt concentration can be read in the figures. In each 
panel, the chain asphericity A at = 6.25 x 10^^ (dashed curve), 7.5 X 10"^ (dotted curve), and 
6.0 X 10^'^ (dash-dotted curve) are plotted to illustrate the chain conformational change. The value 
of A is read from the right y-axis of the plot. 
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effective chain charge of —17.3. This is equivalent to 30.7 monovalent counterions condensing 
on the chain. Our simulations obtained N^-^ = 30.4, very close to the prediction, although 
the chain is flexible and of finite length, which does not follow the assumption of a rigid, 
infinitely long chain. This consistency shows that the theory is a good approximation of 
flexible chains and also that the choice for the condensation threshold = Xb able to 
capture the condensation phenomena. 

In Fig. IHt^b), we observed that A'^^^ remains constant for a much wider range of electric 
field than N^-^. For example, in Cg < 3.75 x 10"^, this number varies very little over 
the studied field strengths, although the electrophoretic mobility (in Fig. [7]) does show some 
variation. Therefore, the decrease of fipe in this salt region is directly related to the stripping- 
off of condensed monovalent counterions from the chain. Moreover, we observed that 
exhibits a two-step decrease at Cg — C*. By referring to the corresponding A curve, we 
can find that the first plateau appears in the weak field region where the chain exhibits 
a unperturbed, globular structure and the condensed tetravalent counterions are wrapped 



within the twisted chain 



28 



29| . The second, small plateau occurs when E is strong and 



the chain is completely stretched with A ~ 1. In this case, the number of the condensed 
tetravalent counterions maintains a constant that is smaller than the globular chain value. 
If the applied field is very strong, such as i? > 1, the tetravalent counterions can be stripped 
off the chain violently; as a consequence, N^^ decreases. The two-step plateaus are smeared 
out when Cg is high. It is because a large number of the tetravalent counterions appears in 
the solution, obscuring the boundary between condensed and non-condensed ions. 

Fig. E^c) shows that in weak fields, there are nearly no colons condensed on the chain when 
Cg < C*. In the high-salt region above the equivalence point C*, the colon condensation 
significantly increases. This is because the effective charge becomes positive, manifested by 
the inversion of chain mobility as seen in Fig. [7](a), and so attracts the colons to the chain. 
After being constant in weak fields, N^^ decreases slightly with E but then displays a large 
peak in strong fields. The peak can be associated with the chain conformational transition 
because tetravalent counterions condensed on a stretched, elongated chain are more exposed 
to the bulk solution than ones condensed/wrapped in a unperturbed, coiled or globular chain 
and hence, attract more colons. The condensed colons constitute an outer layer of the PE 
complex. Therefore, once the chain is fully stretched, the increasing electric field can "blow" 
colons off due to the weak condensation, causing N't^ to decrease drastically. 
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The above results show that the ions do not condense on the chain in a universal way. 
Condensation depends strongly on the ion valency and the chain conformation in electric 
fields. The tetravalent counterions compete with the monovalent counterions. Therefore, 
while A^^4 increases with C^, A^^^ decreases. On the other hand, the colons collaborate 
with the tetravalent counterions; A'^^^^ increases with Cg when Cg > C*, following the trend 
of The non-monotonic behavior of change can be seen in Fig. M^a) where the total 

number of the condensed ions N^^^ is plotted against E at different Cg. Nonetheless, the 




FIG. 9: (a) Total number of condensed ions Nf^^, and (b) total charge of chain, Q{rt = \b), inside 
the condensation region = A^, as a function of E at different Cg- 



total charge inside the condensation region, Q{rt = Xb), (presented in Fig. [9](b)) shows a 
more regular, decreasing behavior. We regarded Q{rt = Xb) as the effective charge and its 
variation compares well with the behavior of the chain mobility fipe shown in Fig. W^a). By 
comparing the results of Fig. |9t^a) and (b), we conclude that it is not the number of ions 
but rather the total charge of the PE-ion complex that is constantly being reduced by the 
strong electric field. 



E. Ion distribution along the chain and the mobilities 

As we have seen and studied, the PE-ion complex is polarized by an external electric field. 
The polarization results from both the migration of the ions condensed on the chain and 
the conformational change. In this subsection, we study how the ion distribution along the 
chain varies with the electric field and see how this distribution is related to the mobilities 
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of the chain and the condensed ions. We associate a condensed ion to the monomer with 
which it is closest, and calculated the mean number n of the condensed ions associated with 
each monomer. The monomers are indexed by i through the chain such that the increment 
of the index follows the sense of the field direction and not necessarily the position along 
the chain. A front end and a rear end of the chain are then defined to be the monomers 
with the smallest and the largest index, respectively. Thus, n{i) is the distribution function 
of the condensed ions. We study four representative cases of Cg- 

The first case considers the salt-free solution. The only ions presented in this case are 
the monovalent counterions. Fig. ITUT a) shows the number distribution of the condensed 
monovalent counterions on the chain, n+i(2). 




monomer ID no. 

(a) 



(b) 



FIG. 10: (a) Distribution of condensed monovalent counterion along chain, n+i(i), at Cg = 0.0 in 
different electric fields E. The value of E is indicated in the figure, (b) Tlie mobility of tlie cfiain 
/ipe, and of the condensed monovalent counterions /i+i, as a function of E. 



In the zero field limit, the n+i curve is fiat in the interior region of the chain and decreases 
near the two chain ends. The decrease can be attributed to the symmetry broken near the 
chain ends because there exists no monomer outside the ends and hence the attractive force 
to condense the counterions is weaker than in the interior. The whole distribution curve is 
symmetric with respect to the middle of the chain. The application of an electric field breaks 
this symmetry by tilting the distribution curve against the field direction, as shown. The 
condensed monovalent counterions accumulate near one end and are depleted near the other, 
leading to a polarization vector pointing in the field direction. The curve shifts downward as 
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E increases. Since the area below the curve represents the number of ions condensed on the 
chain, the downward-shifting shows the decrease of the number of condensed counterions, 
which is consistent with the result in Fig. [Hl^a). The shape of n+i evolves and becomes a 
tangential curve when E is very strong. In order to understand the kinetics of the condensed 
ions, we plot the mobility of the chain /ipg and the mobility of the condensed ions fi'^^ in 
Fig. [TOT b). We found that /i^^ is basically equal to fipe in the weak fields, which shows that 
the condensed counterions reside statically on the chain, /i^^ deviates significantly from /Xpe 
when E goes beyond 0.1. This deviation shows that there is relative motion between the 
condensed monovalent counterions and the chain: The ions glide along the chain. Since 
the chain has finite length, the gliding means that the ions are dynamically, rather than 
statically, condensed on the chain. Once condensed onto the chain, an ion glides along the 
chain backbone and accumulates at the chain end. A condensed ion at the chain end must 
then leave the chain to maintain a steady number of the condensed ions. 

For the second case, consider a salt concentration smaller than the equivalence point, 
Cs = 1.875 X 10"^. According to Fig. [HI there are nearly no colons condensing onto the 
chain. Therefore, we show only the distribution functions for the condensed monovalent 
counterions, n^i{i), and for the condensed tetravalent counterions, n_|_4(z). The results are 
presented in Figs. [TT](a) and (b), respectively. 

Observe that in the zero field limit, the n^i(i) and ^+4(2) curves are both flat in the 
interior and symmetric to the middle of the chain. When an electric field is applied, the 
symmetry is broken, n+i displays a large skewed peak in the interior region of the chain, 
while ?T,_|_4 shows multiple peaks close to the rear end. The number of peaks in n+4 is 3 
for 0.2 < E < 0.5, which corresponds exactly to the number of the tetravalent counterions 
added at this salt concentration, i.e., the tetravalent counterions completely condense. The 
distinct peaks reveals that the ions are localized on the chain, similar to a 1-dimensional 
crystal. This localization happens when appropriate strength of electric field is applied. If 
E is weak, the ions cannot be localized and the ra+4 curve simply tilts. If E is strong, they 
are pushed forcibly to the rear end and only one peak is visible. Similar phenomena have 
been observed for other Cs < C*. The peak of the ra+i curve appears right after the peaks 
of the n_|_4 curve. This non-overlapping of the two species results from mutual exclusion 
between the condensed monovalent counterions and tetravalent counterions due to their 
Coulomb repulsion. Since the external electric force exerted on an ion is proportional to 
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FIG. 11: (a) Distribution of condensed monovalent counterions on chain, n_|_i(i), at Cs = 1.875 x 
10~^(< C*) in different electric fields E. The value of E can be read in the figure, (b) Same as 
(a), but for condensed tetravalent counterions, n+4(z). (c) The mobility of the chain /Xpe, of the 
condensed monovalent counterions /i*^!, and of the condensed tetravalent counterions /Li!^4, as a 
function of E. 



the ion valency, the tetravalent counterions are pushed more strongly by the electric fields 
and move closer to the chain end. The monovalent counterions suffer a weaker force and 
condense just behind the tetravalent counterions. 

We also plot the mobility of the chain and the condensed ions as a function of E in 
Fig. [TTT c) to study the kinetics of the system. We can see that the mobility of the tetravalent 
counterions /i^j_4 coincides with /ipg in small fields and deviates from it when E > 0.5. The 
tight coincidence shows that the tetravalent counterions are tightly bound to the chain and 
move with it. The mobility of the condensed monovalent counterions also takes a value 
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close to /ipe when E < 0.1 but not as close as when Cg = 0.0 (cf. Fig. [101) The tetravalent 
counterions repel the condensed monovalent counterions on the chain, which loosens the 
condensation. Recall that the deviation of fi'^^ from fipe occurs around E = 0.5. This is 
about 4 times the field strength as when /i^^^ begins to deviate from fip^. And the number 
4 is equal to the valency ratio between the two species of the counterions, which can be 
explained by law of friction. Assume that the friction coefficient of the surface is effectively 
The criterion for a condensed ion to glide along the chain surface is governed by Fe > ^Fn, 
where Fe = ZeE, and F„ is the force normal to the chain surface, which is estimated by 
ksTXEZ'^e^ / a"^ . The square of Z appears in F„ because of the matching (or interaction) 
of a Z-valent counterion with Z monomers when it condenses on the chain. Therefore, the 
threshold electric field to glide a condensed counterion is approximately proportional to the 
valency. 

The third case studies the system at the equivalence point C*. At this salt concentration, 
no monovalent counterions or colons condense on the chain; so only the distribution function 
for the tetravalent counterions is shown. The results of n+4(z) and the mobilities for the 
chain and the ions are presented as a function of E in Fig. [T2l 




monomer ID no. E 

(a) (b) 

FIG. 12: (a) Distribution of condensed tetravalent counterions on chain at Cg = 7.5 x 10~^(= CI) 
in different electric fields E. The value of E can be read in the figure, (b) The mobility of the 
chain /ipe and of the condensed tetravalent counterions function of E. 

The behavior of nj^i{i) looks similar to ?t-+i(«) in Fig. [TUT a). The curve tilts, shifts down- 
ward, and eventually shows a tangential distribution as E increases. Since the condensed 
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ions have large valency, the condensation is strong. The chain is charge-neutralized by the 
condensed ions and Fig. [T2r b) shows that fipe is zero when E < 0.5. Above E > 0.5, the 
electric force exerted on the condensed ions can overcome the friction force between the ions 
and the chain surface; hence there is relative motion between the ions and the chain. 

In the fourth case, we study a salt concentration higher than the equivalence point. 
Figs. [TST a) and [TST b) present, respectively, the distributions of the condensed tetravalent 
counterions n+4(z) and of the colons n^i{i) on the chain at Cg = 6.0 x 10~^. The distribution 
of the monovalent counterion n+i{i) is not shown because very few of them condense onto 
the chain. 
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FIG. 13: (a) Distribution of condensed tetravalent counterions on chain, n_|_4(i), at Cs = 6.0 x 
10~^(> C*) in different electric fields E. The value of E can be read in the figure, (b) Same as (a), 
but for condensed coions n_i(i). (c) The mobility of the chain fipe, of the condensed tetravalent 
counterions /U^^4, and of the condensed coions fi'Li, as a function of E. 
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The distribution n+4(i) is similar to the one at C* but the curve is shifted upward by a 
value of about 0.5. This upward-shifting shows that the number of the condensed tetravalent 
counterions exceeds the number needed to neutralize the chain. The chain is overcharged, 
and hence attracts colons onto it. The distribution tilts when E is applied. The tilting 

is not toward the front end but rather toward the rear end and looks similar to n_|„4(i). 
It suggests that the condensation of the colons are, in fact, mediated by the condensed 
tetravalent counterions. The condensation of colon is weak because the mobility for the 
condensed colons is not zero even for the small set E fields simulated, as plotted in 
Fig. [T51(c). The negative value of tells us that the condensed colons migrate constantly 
toward the chain front end. Surprisingly, the number of colons is less elevated at this end. 
This suggests that the chemical potential of the colon condensation is higher at this end. 
The tetravalent counterions condense more numerously at the rear of the chain, lowering 
the chemical potential there. The mobility is identical to /i^g when E < 0.5, which once 
again shows the strong condensation of tetravalent counterions. The tetravalent ions reside 
on the chain and move along with it. The positive value of the mobilities demonstrates 
that the effective chain charge becomes positive due to the over-condensation of tetravalent 
counterions. When E is strong enough to overcome the friction on the chain surface, the 
kinetics for the condensed ions and the chain are decoupled. /i^4 increases but /i^g decreases. 
The magnitudes of the mobilities both increase, /i^^ asymptotically approaches -1, the value 
expected for a colon drifting in the bulk solution. 

We remark that the distributions of condensed ions studied here might be understood 
in the framework of sedimentation. In a dilute solution, sedimentation theory describes a 
barometric density profile p{z) = pQexp{—z/L) for neutral particles, where z is the position 
in the gravitational direction and L = kBT/{mg) is the characteristic length (also equal to 
the mean height of sediments) with m being the particle buoyant mass and g the gravitational 
acceleration. For charged particles, the mean height is shown to extend to ZL owing to 



electrostatic interactions 
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5l| . In our work, the condensed ions are under the action of 



an electric field, which can be analogously regarded as sedimenting in a gravitational field 
with the force mg effectively replaced by the net force ZeE — ^F^ acting on a condensed 
ion. This analogy makes sense only when E is strong enough to displace the condensed 
ions and totally unfold the chain, so that the index i is linearly mapped to z. Following the 
index number in a reverse order, we observe that n^i{i) in Fig. [TU] displays as an exponential 
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function for E > 0.5, except near the two chain ends where some edge effect appears. The 
characteristic length L of the exponential can be seen to decrease with increasing E, which 
follows the depiction of the sedimentation theory. Similar results are found for n^4{i) in 
Figs. [12] and [12] when E is strong. In the second case studied (in Fig. [TT]) . we have two 
species of counterions condensing on the chain. Zwanikken and van Roij have developed 
a mean-field theory for the sedimentation of multicomponent charged colloids, based upon 



a Poisson-Boltzmann approach 



52| A segregation of layering charged colloids of valency 



Zi was demonstrated, in which the order of the sediments, from the bottom to top, follow 



the increasing order of the product ZjLj [51|, 152|. Since Lj = ksT/ {ZicE — and F„ 
is approximately proportional to Zf, the tetravalent counterions have a larger ZiLi than 
the monovalent ones. The theory thus predicts a profile along the chain which has the 
inverse order of what we obtained in the simulations. This difference arises because some 
mechanism cannot occur in mapping to a sedimentation system, for example, the dynamic 
condensation of counterions on a chain. The condensed monovalent counterions are stripped 
off in an electric field, more easily than the tetravalent ones. Consequently, the tetravalent 
counterions stay near the rear end of chain with longer time and the monovalent ones coming 
from the front are then stuck behind the tetravalent, which produces the observed profile. 



F. Dependence of chain mobility on chain length 

Finally we studied the variation of chain mobility fipe with chain length A^, as the applied 
electric E increases. The results are presented in Fig. [TU panel (a) and (b), at two salt 
concentrations, Cs = 0.0 and C*, respectively, for varying from 12 to 768. We have seen 
in Sec. III-C that the magnitude of fipe largely increases when a chain unfolds. 

Observe that the longer the chain length, the weaker the field strength that is needed to 
unfold the chain. Accompanying the unfolding, \fipe\ increases and the mobility reaches a 
constant value in the field region between E = 0.1 and 1. For a stronger E, the magnitude 
of fipe increases again. It is because the condensed ions can be further stripped off the chain, 
which renders the effective chain charge more negative, and thus, the chain drifts faster. The 
electric field required to change the chain mobility depends on the chain length. It provides 
a unique mechanism to electrophoretically separate PEs by A^. Especially at Cs = C*, 
the chains are initially neutralized by the multivalent counterions. They do not drift in 
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FIG. 14: Chain mobility /ipe and degree of chain unfolding Dg versus E (a) in salt-free solutions 
and (b) in solutions at equivalence point C*. Each curve denotes the results run at one chain 
length N and the value of is indicated in the upper panel. 

weak fields. When an appropriate electric field is applied, the longer chains unfold and gain 
mobility, and as a consequence, drift away and are separated from the shorter chains. 

Netz predicted that the unfolding electric field E* should scale as N"'^"^'^ where u is the 
swelling exponent of chain size. We found in Fig. [H] that the chain unfolding, and also the 
mobility change, are not sharp transitions in electric fields. Consequently, it is not easy to 
determine the onset electric field for chain unfolding, Ejj, with a good accuracy. Therefore, 
to verify the Netz' prediction, we calculated, in addition, the electric field Ejjj, which is 
the inflection point on the D^-curve. The results are plotted in Fig. [15] as a function of N, 
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together with E} and E}j. 
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FIG. 15: Critical electric field E* as a function of chain length N for Cg = 0.0 and Cg = Cg. E* 
is determined by the three ways: (1) the polarization energy equal to ksT, (2) the onset point of 
Dg, and (3) the inflection point of Dg. The three estimators are denoted by E"^, Ejj, and Ej^^, 
respectively. 



The three data sets at Cg = 0.0 he on straight hnes in the log-log plot: The critical 
electric field does indeed follow a power-law function. Linear least-square fit yields that 
E*j scales as A^-i-02(4)^ e}^ as A^-o-63(4)^ and E}jj as A^-o-62(3). The scaling law for E} is 



a refinement of our previous work 



19| . It behaves differently to Ejj, which is calculated 



directly from the chain size variation, and thus fails to predict the threshold field to unfold 
a chain. The failure arises from oversimplified setting of the critical polarization energy for 
chain unfolding W*^i to the thermal energy ksT, which should be also a function of chain 
length and salt concentration, as shown in Fig. |6l The two fitting lines for Ejj and Ejjj run 
parallel to each other. Consequently, the ratio of E}jj^ to Ejj is a constant, suggesting that 
either measure of E* is acceptable. At Cg = C*, Ej scales as N~'^'^^^^\ a wrong prediction 
for the unfolding field. For Ejj and Ejjj, markedly different behavior is observed. Only 
the Ejjj data follow a scaling law with the exponent equal to —0.76(2). The onset field 
point Ejj exhibits a crossover from an exponent close to —0.63 at small polymerizations to 
—0.76 when > 100 which agrees with the E}j value. This crossover refiects the fact that 
the transition of chain unfolding becomes sharper when chain length increases, as seen in 
Fig.[Hb). 



30 



To check Netz' theory [3, Q] , we calculate the swelling exponent v of chains in the zero 
field limit through the scaling relation i?^ ~ N'^^ . Our simulations show that v = 0.92(1) at 
Cs = 0.0 and u = 0.33(1) at Cg = C*. Therefore, according to Netz, the critical field should 
scale as A^~^-2^(2) and A^"'^-^''^^) , respectively. However, the scaling law obtained in Fig. [TS] 
differs from the prediction. The modified Netz' theory that considers PEs as ellipsoidal 
objects predicts that the critical electric field scales as V^"^/^ where V is calculated from 



the eigenvalues of the gyration tensor 19(]. Our refining results show that V ~ A^2.44(i) 
Cs = 0.0 and V ~ iV^'^'^*^^) for Cs = C* in zero fields. The modified Netz' theory hence 
predicts E* ~ iV-i-22(i) and E* ~ iV-o-56{i)^ which are, again, inconsistent with the results 
of FigJT^ The differences are so important that a new and detailed understanding of the 
mechanism of chain unfolding is necessary. It is definitely a topic worthy to be investigated 
in the future. 

We remark that Netz' original work does not verify the scaling of E* with chain length. 
Moreover, the PE was collapsed by monovalent counterions (g = 1) due to the assumption 
of a large Coulomb coupling constant S = Xsq^ /u-, ranging between 5 and 30. The situation 
is equivalent to a collapsed PE in a (l:l)-salt solution at very low temperature because the 
Bjerrum length is large, with a value lying between 5cr to 30o". In the derivation, no 
multivalent counterions were considered, neither was the competition between multivalent 
ions and monovalent ones as happens in reality. Therefore, the disagreement with our results 
arises from these simplified assumptions in the model. In the theory, the critical polarization 
energy is equal to ksT, which has to be improved. Moreover, the transition is assumed to 
be sharp. Our simulations reveal a more complicated story: The chain conformational 
transition is a continuous change at Cs = 0.0, while it is a more discontinuous transition 
at Cs = C*, particularly when is large. Therefore, the onset critical field E^ and the 
inflection critical field Ejjj follow different scaling laws in different solutions. The ratio 
between the two estimators describes the sharpness of the transition and theory should 
account for its deviation from unity. 

In our previous work 20(], the critical electric field in trivalent salt solutions at the equiv- 
alence point was studied through the infiection-point method. The scaling exponent was 
found to be —0.77(1), which is consistent with what we obtain here. This seems to suggest 
that the infiection critical field for trivalent and for tetravalent salts follows a similar scaling 
law. Liu et al. have recently studied the scaling law of DC unfolding fields in monovalent. 
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divalent, and trivalent salt solutions by simulations [2l|. The exponents reported depended 
on the salt valency. The result for trivalent salt (—0.64) at equivalence point is not consistent 
with ours. Since the chain length is short in their study (only up to = 192 in comparison 
with ours = 768) and the inflection point is ambiguous (see Fig. 3 in the paper), their 
results are not reliable. However, the study gives rise a relevant and interesting question: 
Does the critical electric field follow a different scaling law when the valency of salt is small, 
such as monovalent and divalent salt, as claimed? It necessitates a detailed and precise in- 
vestigation in the future. According to our work, the unfolding electric field shares a similar 
scaling exponent in trivalent and tetravalent salt systems but the prefactor in the scaling 
law is different. Many phenomena in both of the systems occur in a similar way, such as ion 
distributions and mobility changes. 

We study the tetravalent system because large ion valency can give a stronger response 
to applied electric fields, which makes observing the effects easier. From the point of view 
of applications, it is very important to understand the role of salt valency in PE solutions in 
electric fields. The selection of salt valency depends on the research context and goal. We 
hope that the information reported here can be helpful in the development of techniques for 
molecular separation and in the design of functionalized micro/nano-fluidic devices. 



IV. CONCLUSIONS 



In this work, we performed molecular dynamics simulations to investigate the conforma- 
tional and electrophoretic properties of chains in tetravalent salt solutions subject to electric 
fields. Our results show that chain size depends strongly on the salt concentration Cg and 
that under the action of electric fields, the chain shape can be altered. When the field 
strength is stronger than a critical value E* , the chains are largely extended to an elongated 
structure. Two estimators and E*jj were used to calculate E* through equating Wpoi to 
ksT and identifying the onset point of Dg, respectively. The obtained values show that E* 
is a non-monotonic function of C^, and the maximum value appears at the equivalence point 

= C*. The dipole moment shows that chain polarization displays a linear dependence on 
the electric field up to E}^. 

The salt concentration has a strong influence on the electrophoretic mobility of the chain 
and ion distributions. In weak electric fields, the chain mobility /Xpe is negative for < C*, 
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whereas it is positive when Cs > C*. The latter demonstrates the sign inversion of the 
effective chain charge. The mobihty of the condensed tetravalent counterions /i^4 is identical 
to fipe because the condensed ions are tightly bound onto the chains and move together in 
the electric fields. When the applied electric field E is stronger than Ejj, the chains unfold 
such that /i^j_4 and fipe are no longer identical. In an even higher E field, becomes 
positive and fipe becomes negative: There is a relative motion between the chains and the 
condensed tetravalent ions. The number of the condensed ions shows that part of the 
condensed tetravalent ions are stripped-off from the chains due to the strong electric fields. 
The behavior of the effective chain charge Q is consistent with the chain mobility found. 

Moreover, we studied in detail the distribution of different species of ions condensed on 
the chains. The counterions are dragged toward the rear of the chain, due to the polarization. 
For Cs < C*, the condensation profile displays that the tetravalent counterions condense at 
the chain rear, followed by the condensed monovalent ones right behind the tetravalent, close 
to the chain center, owing to the strong electrostatic repulsion between the counterions. For 
Cs ^ C*, monovalent counterions are totally expelled from the chain, which leaves only the 
condensed tetravalent counterions, which overcharge the surface of the chains due to the 
excessive condensation. Consequently, the colons are attracted to the chains and condensed. 
This process is mediated by the condensed tetravalent counterions because a similar profile of 
colon distribution to the tetravalent one is found. An analogy of our system to sedimentation 
problems was used to explain the ion condensation profile. 

Finally, we investigated the dependence of chain mobility and unfolding transition on the 
chain length N. The estimator Ejjj was considered through the calculation of the inflection 
point of Dg. At Cs = 0.0, the transition of chain size is not sharp. But Ejj and Ejjj follow 
a similar scaling law A^~o-^3(^) with different prefactors. However, at the equivalence point, 
the unfolding transition becomes sharper. Ejjj scales as A^~*'-'^^(2)_ crossover was observed 
in Ejj, which converges asymptotically to Ejjj when is large. The scaling laws obtained 
here are significantly different to both Netz' prediction and the modified theory. Noticeably, 
Ej fails to predict the unfolding field, which shows that the critical polarization energy is 
not simply /c^T. It hence necessitates a further investigations to explore the mechanism of 
chain unfolding in electric fields and the associated mobility and electrokinetics changes. 
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